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Abstract 

This paper presents a theoretical study of blood flow through a tapered and overlapping 
stenosed artery under the action of an externally applied magnetic field. The fluid (blood) 
medium is assumed to be porous in nature. The variable viscosity of blood depending on 
hematocrit (percentage volume of erythrocytes) is taken into account in order to improve 
resemblance to the real situation. The governing equation for laminar, incompressible 
and Newtonian fluid subject to the boundary conditions is solved by using a well known 
Frobenius method. The analytical expressions for velocity component, volumetric flow rate, 
wall shear stress and pressure gradient are obtained. The numerical values are extracted 
from these analytical expressions and are presented graphically. It is observed that the 
influence of hematocrit, magnetic field and the shape of artery have important impact on 
the velocity profile, pressure gradient and wall shear stress. Moreover, the effect of primary 
stenosis on the secondary one has been significantly observed. 

Keywords: Overlapping Stenoses, Tapered Artery, MHD Flow, Porous Vessel, Hemat- 
ocrit, Frobenius Method 



1 Introduction 

Many cardiovascular diseases such as due to the arterial occlusion is one of the leading cause 
of death world wide. The partial occlusion of the arteries due to stenotic obstruction not 
only restrict the regular blood flow but also characterizes the hardening and thickening of the 
arterial wall. However, the main cause of the formation of stenosis is still unknown but it is well 
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established that the fluid dynamical factors play an important role as to further development 
of stenosis. Therefore, during the past few decay several studies were conducted by Young [1], 
Young and Tsai [2, 3] to understand the effects of stenosis on blood flow through arteries. 



Nomenclature 




z 


Axial distance 


r 


Radial distance 


R(z) 


Arterial wall 


d 


Distance from inlet of onset of stenosis 


u 


Axial velocity component 


u 


Non-dimensional axial velocity 


H{r) 


Coefficient of viscosity of blood at a radial distance r 


Ho 


Viscosity coefficient for plasma 


h(r) 


Hematocrit at a distance r 


H 


Maximum hematocrit at the center of the arterial segment 


m(>2) 


Shape parameter of hematocrit 


Rq 


Radius of the pericardial surface of the normal portion of the arterial segment 


I 


Length between throat of two stenosis stenoses 


L 


Length of arterial segment 


Bo 


Magnetic field strength 


a 


Electrical conductivity 


k 


Porous permeability parameter 


M 


Hartmann number 


V 


Blood pressure 


dp 

dz 


Non-dimensional axial pressure gradient 


Q 


Volumetric flow rate 


Q 


Non-dimensional volumetric flow rate 




Wall Shear stress 


f 


Non-dimensional wall shear stress 



Tu and Deville [4] investigated pulsatile flow of blood in stenosed arteries. Misra and Shit 
[5, 6] studied in two different situations on the blood flow through arterial stenosis by treating 
blood as a non-Newtonian ( Herschel-Bulkley ) fluid model. It is generally well known that 
blood, being a suspension of red cells in plasma, behaves like a non-Newtonian fluid at low 
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shear rates [7, 8]. However, Misra and Chakravarty [9] and Shit and Roy [10] put forwarded a 
mathematical analysis for the unsteady flow of blood through arteries having stenosis, in which 
blood was treated as a Newtonian viscous incompressible fluid. 

The hemodynamics associated with a single stenotic lesion are significantly affected by the 
presence of a second lesion. In many situations there are evidences of the occurrence of the 
multiple or overlapping stenosis such as the patients of angiograms. Misra et al. [11] conducted 
a theoretical study for the effects of multiple stenosis. An experimental study of blood flow 
through an arterial segment having multiple stenoses were made by Talukder et al. [12]. The 
effects of overlapping stenosis through an arterial stenosis have been successfully carried out 
analytically as well as numerically by Chakravarty and Mandal [13] and Layek et al. [14] 
respectively. However, all these studies are restricted in the consideration of magnetic field and 
the porous medium. 

Since blood is electrically conducting fluid, its flow characteristics is influenced by the ap- 
plication of magnetic field. If a magnetic field is applied to a moving and electrically conducting 
fluid, it will induce electric as well as magnetic fields. The interaction of these fields produces 
a body force per unit volume known as Lorentz force, which has significant impact on the flow 
characteristics of blood. Such an analysis may be useful for the reduction of blood flow during 
surgery and Magnetic Resonance Imaging ( MRI ). The effect of magnetic field on blood flow has 
been analyzed theoretically and experimentally by many investigators [15] - [18] under different 
situations. Shit and his co-investigators [19] - [23] explored variety of flow behaviour of blood 
in arteries by treating Newtonian/ non- Newtonian model in the presence of a uniform magnetic 
field. 

Very recently, the studies of blood flow through porous medium has gained considerable 
attention to the medical practitioners/ clinicians because of its enormous changes in the flow 
characteristics. The capillary endothelium is, in turn, covered by a thin layer lining the alveoli, 
which has been treated as a porous medium. Dash et al. [24] considered the Brinkman equation 
to model the blood flow when there is an accumulation of fatty plaques in the lumen of an 
arterial segment and artery-clogging takes place by blood clots. They considered the clogged 
region as a porous medium. Bhargava et al. [25] studied the transport of pharmaceutical species 
in laminar, homogeneous, incompressible, magneto-hydrodynamic, pulsating flow through two- 
dimensional channel with porous walls containing non-Darcian porous materials. Misra et al. 
[20, 26] presented a mathematical model as well as numerical model for studying blood flow 
through a porous vessel under the action of magnetic field, in which the viscosity varies in the 
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radial direction. 

Hematocrit is the most important determinant of whole blood viscosity. Therefore, blood 
viscosity and vascular resistance ( due to the presence of stenosis ) affect total peripheral re- 
sistance to blood flow, which is abnormally high in the primary stage of hypertension. Again 
hematocrit is a blood test that measures the percentage of red blood cells present in the whole 
blood of the body. The percentage of red cells in adult human body is approximately 40 - 45 
% [27]. Red cells may affect the viscosity of whole blood and thus the velocity distribution 
depends on the hematocrit. So blood can not be considered as homogeneous fluid [14]. Due to 
the high shear rate near the arterial wall, the viscosity of blood is low and the concentration of 
cells is high in the central core region. Therefore, blood may be treated as Newtonian fluid with 
variable viscosity particularly in the case of large blood vessels. 

The present study is motivated towards a theoretical investigation of blood flow through a 
tapered and overlapping stenosed artery in the presence of magnetic field. The study pertains 
to a situation in which the variable viscosity of blood depending upon hematocrit is taken into 
consideration. The present model is designed in such a way that it could be applicable to both 
converging/ diverging artery depending on the choice of tapering angle a. Thus, the study will 
answers the question of mechanism of further deposition of plaque under various aspects. 



2 Mathematical modelling of the problem 

We consider the laminar, incompressible and Newtonian flow of blood through axisymmetric 
two-dimensional tapered and overlapping stenosed artery. Any material point in the fluid is 
representing by the cylindrical polar coordinate (r, 9, z), where z measures along the axis of the 
artery and that of r and 9 measure along the radial and circumferential directions respectively. 
The mathematical expression that corresponds to the geometry of the present problem is given 
by 



R{z) = R 



1.0 



11.0 
323 



47.0 



= Roa(z) 



f{z -<£) + j^l 2 {z - df - l(z - df + (1.0/3.0)(2 - df a{z) 
48.0 J 

3/ 

where d < z < d + — 
elsewhere 



(1) 



where the onset of the stenosis is located at a distance d from the inlet, W the length of the 



J > 2 



overlapping stenosis and I representing the distance between two critical height of the stenoses. 
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The expression for a(z) is responsible for the artery to be converging or diverging depending on 
tapering angle a has the form 

a(z) = 1 + z tan(a). (2) 




1 2 3 4 5 



Fig. 1 Schematic diagram of the model geometry. 



We assumed that blood is incompressible, suspension of erythrocytes in plasma and has uniform 
dense throughout but the viscosity /i(r) varies in the radial direction. According to Einstein's 
formula for the variable viscosity of blood taken to be 

/i(r)=/i [l + 0h(r)], (3) 

where /xq is the coefficient of viscosity of plasma, /3 is a constant (whose value for blood is 
equal to 2.5) and h(r) stands for the hematocrit. The analysis will be carried out by using the 
following empirical formula for hematocrit given by in Lih [29] 

Mr)=#[i-(fn (4) 



in which R represents the radius of a normal arterial segment, H is the maximum hematocrit 
at the center of the artery and m(> 2) a parameter that determines the exact shape of the 
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velocity profile for blood. The shape of the profile given by Eq.(4) is valid only for very dilute 
suspensions of erythrocytes, which are considered to be of spherical shape. 

According to our considerations, the equation that governs the flow of blood under the 
action of an external magnetic field through porous medium may be put as 

dz r dr +aB » u+ k U ~^ {b) 

where u denotes the (axial) velocity component of blood, p the blood pressure, a the electrical 
conductivity, k the permeability of the porous medium and -Bo is the applied magnetic field 
strength. 

To solve our problem, we use no-slip boundary condition at the arterial wall, that is, 

u = at r = R(z). (6) 

Further we consider axi-symmetric boundary condition of axial velocity at the mid line of the 
artery as 

du 

— = at r = 7) 
or 

3 Analytical solution 

In order to simplify our problem, let us introduce the following transformation 

r 

£ = TT (8) 

With the use of the transform defined in (8) and the equations (3) and (4), the governing 
equation (5) reduces to 

with ai = 1 + a 2 , a 2 = /3H, k = -j^ and M 2 = a{B °* o)2 ■ 
Similarly the boundary conditions transformed into 

u = at (10) 

du 

and — = at f = (11) 



The equation (9) can be solved subjected to the boundary conditions (10) and (11) using Frobe- 
nius method. For this, of course, u has to be bounded at £ = 0. Then only admissible series 
solution of the equation (9) will exists and can put in the form 

oo I?2^P oo 

« = #E4f + p* L E s *r +2 , (12) 

i=0 4 °i/ i o i=0 

where K, A { and Bi are arbitrary constants. 

To find the arbitrary constant K, we use the no-slip boundary condition (10) and obtained as 

K = -p±*% — . (13) 

Substituting the value of u from equation (12) into equation (9) and we get 

oo oo oo 

l)(oi - a 2 r)Af" 2 + E<(ai - (m - ljaD^- 2 - (M 2 + y) E^ 

i=0 i=0 K i=0 

oo f?2(ip oo oo 

+ El^ i+ 1 + T^fEG + + 2 )(«i - ^OBiC + E(< + 2)(ai " (™ " l)«2D^r 

i=0 fc 4a l/^0 L i= o i=0 



oo oo 



p2 dp 

-(m 2 + ^)E ^e +2 + E ^ +m+2 l = t^( 14 ) 



^ i=0 i=o k 



4ai/i 



Equating the coefficients of if and other part in equation (14) we have, 



^ - i)( 0l - a 2 r)Af- 2 + e*(«i - (™ - iKr)Ar 2 - (M 2 + ^ ) e 

i=0 i=0 ft i=0 

oo 

+E^e +m = o (is) 

i=0 ft 

and 

oo oo 

E(* + i)(i + 2)( ai - a 2 r)5,e + E(* + 2 )(«i - (™ - i^rw 

i=0 i=0 

-(m 2 + ^)E B t e +2 + e ^e +m+2 = i (i6) 

ft i=0 i=0 ft 

Hence the constants and B, L are obtained by equating the coefficients of and from 
both side of equations (15) and (16) respectively and can be put in the form 

q 2 (l + i) 2 A l+1 ^ m - ma 2 A l+1 _ m + (M 2 + f)A-i - f A-i-m 



ai(l + i) 



2 



(17) 



a 2 (3 + *) 2 - a 2 mB i+1 _ m + (M 2 + f )B i _ 1 - ^B,_ x _ m 

ai(3 + ?) 2 



with 

A = B = 1. 

Substituting the expression for K in the equation (12), we have 



(19) 



u 



e ^(i) i+2 e - e ^e +2 e ^(f ) s 

i=0 i=0 i=0 i=0 



4ai// 

The average velocity m has the form 



EA(£)* 



(20) 



«o = 



8/i V^/o 



(21) 



where is the pressure gradient of the flow field in the normal artery in the absence of 

magnetic field. The non-dimensional expression for u is given by 



u 



u 



dp 
dz 



oo oo oo oo 

E ^(f)* +2 E M* - E B£+ 2 E A^y 

i=0 m i=0 i=0 i=0 "° 



u o «i (So 



(22) 



i=0 



The volumetric flow rate across the arterial segment is given by 

R 
R 

Q= J 2irRotu(Odt (23) 
o 

Substituting u from Eq.(20) into Eq.(23) and then integrating with respect to £, we obtain 



Q = 



tu ,e n (<+2) 



fe) I+2 - Bi (£) i+4 CO 



+4) ^AUJ 



i=0 



i=0 



i=0 



2ai/i 



EA(f) 2 



(24) 



i=0 



If Qo be the volumetric flow rate in the normal portion of the artery, in the absence of magnetic 
field and porosity effect, then 

_ ttR 3 dp 



(25) 



Therefore, the non-dimensional volumetric flow rate has the following form 



O 4^ 
Q=¥-- dz 



oo A-(S-) i+2 oo s-(_R ) i+4 oo 



A<- £ M.r ,I A (i +2 ) 

«=0 «=0 



E 



i=0 i=0 



Q° ai(f)o 



EM0 



(26) 
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If the flow is steady and no outward/inward flow takes place through the arterial segment, then 
the mass flux is constant and hence Q — 1. The expression for pressure gradient from (26) can 
be put as 



dp 
dz 



,dp _ 

^ (So 4 



oo „ oo Ai(^-) i+2 oo B4(S-) i+A oo 

^ n n ^R ) Ac (i+2) A (i+4) Ac^MflJ 

t=0 1=0 t=0 1=0 



(27) 



The wall shear stress on the endothelial surface is given by 

du 



- fi(r) 



dr 



(28) 



J r=R(z) 



Substituting u from Eq.(20) into Eq.(28), we obtain 



tr 



^o[l + Wl-(i) m } 



g r 2 g <A(S) i_1 - g (< + mm i+1 g A(f ) 



i=0 



i=0 



i=0 



Ro< 



-Ro- 



4ai 



(29) 

If rjv = — "^(^f)o be the shear stress at the normal portion of the arterial wall, in the absence 
of magnetic field, the non-dimensional form of the wall shear stress is given by 



tr 
tn 
1 dp 



R 



dz 



, )[i +pH(i -(— m- „ 

2ai i? EA(fyU=-o 



-E( < + 2 )^ 



*=o '*<>' 



p oo p 

E^(#r +2 E^(#) 

-Kq i=0 -Ko 



i-1 



i=0 







Kq 



i=0 



(30) 



4 Results and Discussion 

In the previous section we have obtained analytical expressions for different flow characteristics 
of blood through porous medium under the action of an external magnetic field. In this section 
we are to discuss the flow characteristics graphically with the use of following valid numerical 
data which is applicable to blood. To continue this section we have use the following standard 
values of physical parameter: 

/ = 2.0, d = 0.5, a = 0.09, H = 0.2, k = 0.25, (3 = 2.5, L = 5, M — 2.5, 
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Fig. 2 indicates the different locations of the stenosis in the axial direction. In Fig. 2, z — 0.5 
and z = 3.5 correspond to the onset and outset of the stenosis and z = 1 and z = 3 represent 
the throat of the primary stenosis and secondary stenosis respectively. It is interesting to note 
from this figure that z = 2 is the location where further deposition takes place and hence it is 
known as overlapping stenosis. 
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Fig. 2 Locations of the points z=0.5, 1.0, 2.0, 3.0, 3.5. 



The variation of axial velocity at different axial position along the radial direction are shown 
in Fig. 3. We observe from this figure that the velocity is maximum at the central line of the 
vessel for all position of z. Among all these positions, the velocity is high at the throat of the 
primary stenosis and low at the onset of the overlapping stenosis. However, the central line 
velocity at the throat of secondary stenosis is about 30 % less than the primary stenosis. But 
the central line velocity at z = 2 (between the throat of two stenosis) suddenly falls about 55 
% than that of secondary one. This observation may leads to the flow circulation zone and may 
causes further deposition of plaque. Fig. 4 depicts the variation of axial velocity for different 
length between the throat of two stenoses. It has been observed that the velocity at the throat of 
the stenoses significantly increases with the increase of I. Thus, the effect of the shape of stenosis 
has important role on the flow characteristics. Similar is the observation from Fig. 5 that the 
axial velocity decreases as the tapering angle a increases. Fig. 6 illustrates the variation of axial 
velocity at the throat of the secondary stenosis for different values of the Hartmann Number M. 
We observe that the axial velocity significantly decreases at the central line of the artery with 
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the increase of the magnetic field strength. While the velocity in the vicinity of the arterial wall 
increases with the increasing values of M in order to maintain constant volumetric flow rate. 
It is also well known that when a magnetic field is applied in an electrically conducting fluid 
(here for blood) there arises Lorentz force, which has a tendency to slow down the motion of 
the fluid. It has been observed from Fig. 7 that the axial velocity near the central line of the 
channel increases with the increase of the permeability parameter k, while the trend is reversed 
in the vicinity of the arterial wall. This phenomena is noticed because of the permeability 
parameter k is depend as the reciprocal of the permeability of the porous medium. Fig. 8 gives 
the distribution of axial velocity for different values of the hematocrit H. We note from Fig. 8 
that the axial velocity decreases at the core region of the artery with the increase of hematocrit 
level H, whereas the opposite trend is observed in the peripheral region. This fact lies with in 
the hematocrit as the blood viscosity is high in the core region due to the aggregation of blood 
cells rather than low viscosity in the plasma near the arterial wall. 

Figs. 9-12 illustrate the variation of pressure gradient ^ along the length of the stenosis 
for different values of the physical parameters of interest. Fig. 9 shows that the axial pressure 
gradient increases with the increase of magnetic field strength. We have already observed that 
the Lorentz force has reducing effect of blood velocity, so as more pressure is needed to pass 
the same amount of fluid under the action of an external magnetic field. However, the opposite 
trend is observed in the case of porous permeability parameter k as shown in Fig. 10. It has 
been seen from Fig. 11 that the pressure gradient increases with the increase of the hematocrit 
H . It indicates from this figure that when the aggregation of blood cells increase at the core 
region that is hematocrit H is high, more pressure gradient is needed to pass the same amount 
of the fluid through the stenotic region. It is interesting to note from these three figures that 
the magnitude of the pressure gradient is high enough at the throat of the primary stenosis 
than that of the secondary one. But from Fig. 12, we observed that at the throat of secondary 
stenosis, the pressure gradient is high in comparison to the throat of the primary stenosis. This 
happens due to the increasing of the tapering angle a. Therefore, in the case of diverging artery 
more pressure is needed as the flow advances in the downstream direction. 

Figs. 13 and 14 give the distribution of the wall shear stress for different values of the 
hematocrit H and tapering angle a. We observe from Fig. 13 that the wall shear stress increases 
as the hematocrit H increases. One can note from this figure that the wall shear stress is low at 
the throat of the secondary stenosis as well as at the downstream of the artery. It is generally 
well known that at the low shear stress region mass transportation takes place and thereby 



11 



occurs further deposition. However, it is interesting to note from Fig. 14 that the wall shear 
stress decreases significantly with the increasing values of the tapering angle a. Moreover, the 
magnitude of the wall shear stress is same at both the throat of the stenosis in the absence of 
tapering angle a. Therefore, we may conclude that there is a chance of further deposition at 
the downstream of the diverging artery. 

5 Conclusions 

A theoretical study of blood flow through overlapping stenosis in the presence of magnetic field 
has been carried out. In this study the variable viscosity of blood depending on hematocrit and 
the has been treated as the porous medium. The problem is solved analytically by using Frobe- 
nius method. The effects of various key parameters including the tapering angle a, percentage 
of hematocrit H, the magnetic field and permeability parameter k are examined. The main 
findings of the present study may be listed as follows: 

• The effect of primary stenosis on the secondary one is significant in case of diverging artery 
(a > 0) 

• The flow velocity at the central region decreases gradually with the increase of magnetic 
field strength. 

• The permeability parameter k has an enhancing effect on the flow characteristics of blood. 

• At the core region, the axial velocity decreases with the increase of the percentage of 
hematocrit H. 

• The hematocrit and the blood pressure has a linear relationship as reported in [28]. 

• The lower range of hematocrit may leads to the further deposition of cholesterol at the 
endothelium of the vascular wall. 

Finally we can conclude that further potential improvement of the model are anticipated. Since 
the hematocrit positively affects blood pressure, further study should examine the other factors 
such as diet, tobacco, smoking, overweight etc. from a cardiovascular point of view. More over 
on the basis of the present results, it can be concluded that the flow of blood and pressure can 
be controlled by the application of an external magnetic field. 
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Fig. 3 Velocity distribution in the radial direction at different axial position z, when H 

M = 2.5, (3 = 2.5, a = 0.09 and k = 0.25. 




Fig. 4 Variation of axial velocity along axial direction for different values of /, when H 

M = 2.5, a = 0.09 and k = 0.25. 
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Fig. 5 Variation of axial velocity at z — 2.0 for different values of a, when H = 0.2, M 

/3 = 2.5, and k = 0.25. 
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Fig. 6 Velocity distribution at z — 2.0 for different values of M, when H = 0.2, (3 = 2. 

a = 0.09 and k = 0.25. 
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Fig. 7 Variation of axial velocity for different values of the permeability parameter k at 
z = 2.0, when H = 0.2, M = 2.5, (3 = 2.5 and a = 0.09. 
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Fig. 8 Variation of axial velocity in the radial direction for different values of H z = 2.0, when 

M = 2.5, (3 = 2.5, a = 0.09 and k = 0.25. 



18 



dp 

dz 




z 
L 



;. 9 Variation of pressure gradient for different values of M when H = 0.2, f3 

a = 0.09 and k = 0.25. 
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g. fO Variation of pressure gradient (-£) with z for different values of when if 

M = 2.5, = 2.5, when a = 0.09. 
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11 Variation of pressure gradient (^) for different values of H when M = 2.5, (3 



> dz 

a = 0.09 and k = 0.25. 
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g. 12 Variation of pressure gradient (^) with z for different values of a when H 

M = 2.5, (3 = 2.5, and k = 0.25. 
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Fig. 13 Distribution of wall shear stress r along with z for different values of H when k = 0.25, 

M = 2.5, (3 = 2.5, and a = 0.09 




Fig. 14 Distribution of wall shear stress r along with z for different values of a when H = 0.2, 

(3 = 2.5, M = 2.5 and k = 0.25. 
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